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Abstract. The first 52 fs of a time evolution of the electron density in OCS after an interaction with an 
intense sub 10 fs laser pulse are studied using the time-dependent density functional theory. The nuclear 
motion in this linear trimer is simulated by the classical molecular dynamics method. Laser fields of 
intensity 10^^ W/cm^ and 10^^ W/cm^ are used. Details of the laser induced changes of the structure, as 
well as the ionization rate are sensitive to the applied field intensity and its polarization. It is found that 
under suitable conditions the OCS molecule bends soon after an interaction with a laser pulse. A deviation 
from the linear geometry of up to 23.6° and charged ions of up to +3 are observed. The time evolution of 
electric dipole moments and the time-dependent electron localization function (ELF) are also studied. 

PACS. 31.15.ee Time-dependent density functional theory - 33.80.-b Photon interactions with molecules 
- 36.40.Qv Stability and fragmentation of clusters 

1 Introduction intervals below 10"'^'^ s as ultra-short laser pulses. Ex- 
perimental and theoretical studies have been shown that 
An interaction between nanoparticles and strong laser fields the multi-photon ionization mechanism dominates a laser- 
is now an active area of research |ll2l3l4| . Femtosecond molecule interaction at intensities < lO^^ W/cm^ for laser 
pulsed lasers are able to produce high-intensity fields and frequencies in the near infrared part of the spectrum. The 
they are used to initiate and investigate various processes tunnel and multi-electron dissociative ionizations prevail 
in nanoparticles [Tl^ . It is customary to caU intensi- fo'' intensities greater than IQi^ W/cm^. Fields be- 
ties above lO^^ W/cm^ as intense laser fields, and time lO^^ W/cm^ are not intense enough to induce the 
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multiple ionization of neutral molecules. Research efforts ifications of the shape in chemical processes and under 
are now oriented towards attosecond laser pulses and in- strong fields. It has been shown that the OCS molecule 
tensities of 10^'* W/cm^ [31. When atoms, molecules, or at high energies close to dissociation exhibits chaotic dy- 
clusters, accelerated to a high kinetic energy, impinge on namics [T5]. The Coulomb explosion of carbonyl sulphide 
thin foils they sometimes undergo through the process of has been investigated using ion-momentum imaging tech- 
Coulomb explosion. It is also possible to achieve similar niques [H]. The laser pulses of 55 fs, up to 2 x 10^^ W/cm^ 
conditions if a nanoparticle collides with an accelerated of intensity, and wave length of A = 790 nm have been 
beam of other particles, or interacts with intense laser used. The shape of the exploding OCS molecule has been 
fields [516] . In these processes often occurs the multiple determined. Measured ion momentum maps have been 
ionization of a targeted nanoparticle. This nanoparticle compared to those obtained using the Monte Carlo sim- 
then explodes as a result of the strong Coulomb repul- ulation. The ion momentum maps have been shown that 
sion. One or several charge centers may form in an initial 0"^+ and S'^+ exist along the laser polarization direction, 
nanoparticle when the loss of electrons occurs under laser The central C'^"^ exhibits a four-lobe structure which indi- 
light. Products of the Coulomb explosion are high-energy cates that a non-linear geometry of OCS develops in the 
electrons and ions, and under suitable conditions even X- Coulomb explosion. A peak bond angle of 170° and bonds 
rays. A behavior of fragments in the Coulomb explosion stretching have been measured. The 0-C and C-S bonds 
sheds light on interactions in an initial nanoparticle. The doubled in an interaction with laser pulses. In a recent 
Coulomb explosion imaging is successfully used to deter- work the stronger field of ~ 10^^ W/cm^ has been used 
mine an equilibrium structure of nanoparticles. to investigate the multiple ionization of OCS [14] . The 



The carbonyl sulphide is an important material present 
in the Earth's atmosphere, volcanic gases, ice cores in 
Antarctica, and the interstellar medium. It catalyzes the 
formation of biological molecules and it has been sug- 
gested that OCS played an important role in the origin 
of life on our planet [71 . The charged clusters of carbonyl 
sulphide have been also studied [819] . In its ground state 
the OCS molecule is linear, but it is bent in the first ex- 



three dimensional covariance mapping technique has been 
applied. 

The time-dependent density functional theory (TDDFT) 
|15ll6ll7j is a powerful technique for a description of ex- 
citations in nanoparticles. It is possible to use this method 
to study particle dynamics under strong laser pulses [1611811912012112!^ . 
The time-dependent Schrodinger equation is replaced by 
the time-dependent Kohn-Sham one (atomic units are used) 



d 



cited state |10lllj . In addition, the OCS molecule consists Qt 



1 



i77:4(r,i)= -- V +i'A's(r, t) 0»(r, i). (1) 



2 



of three different atoms. Therefore, in comparison with In this equation uj^5(r, i) = (r, -t-Uj.c(r, t) -|-Uea;t(r, t), 
homonuclear trimers, it is more prone to various mod- where vh is the Hartree potential, Vxc is the exchange- 
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correlation potential, and Vf-^t is an external potential, light. In the theoretical analysis of the lascr-molcculc in- 
The particle density n is calculated from teraction it is difficult to treat on the same level situations 

where the particles in the molecule are bonded together 

n(r,i)=^0:(r,i)0.(r,t), (2) 

* and where nuclear densities are delocalizcd. The classical 

where N is the number of occupied states. The laser field ^ ^ , ^ i • i i i i i -i- r i- 

treatment of nuclei brafces when a probability for a dis- 

enters in an external term of the Kohn-Sham potential. . . „ , , , . , , , , . , . , n k i ah i 

sociation of the molecule m the laser held is high |f 7l25j . 

In the classical dipole approximation the potential of the „, , , . r .-i t i . i ; • 

ihe proper treatment of the coupling between electronic 

laser field is 

and nuclear motion under strong laser light should be car- 

JV 

V{v,t) = If{t)sm{ujt)^r ■ a, (3) "ed out using the multi-component TDDFT [26|27j. In 

this theory both nuclei and electrons are treated as quan- 

where / denotes the intensity of the field, oj is the fre- 

tum particles. However, this approach is still in develop- 

quency, and the a is the polarization of the light. The 

ment. The multi-component TDDFT is computationally 

function f{t) determines the shape of the laser pulse. The 

time consuming and it has been performed only for small 

laser potential is a strong term and it is necessary to solve 

systems. 

the Kohn-Sham equations in a non-perturbative regime. 

In this work; mixed techniques of the quantum me- 

The Born-Oppenhcimer approximation is often applied 

chanical time-dependent density functional theory and the 

in the analysis of the molecules. It assumes that the mo- 
classical dynamics method are used to study the heteronu- 

tion of nuclei and electrons can be separated. Nuclei are 

clear linear trimer OCS after an irradiation by an intense 

much heavier and therefore their motion often can be ne- 

and short laser pulse. The results show a bending of the 

glected when electrons move. However, in some physical 

molecule and an elongation of its bonds, as well as ion- 
situations the Born-Oppenhcimer approximation breaks, 

izations, intensive energy exchanges and charge redistri- 

and it is important to treat the coupled dynamics of elec- 

butions. In the following the computational method is de- 

trons and nuclei. One example of a such coupled motion 

scribed in Sec. II. Results and discussion are presented in 

occurs for molecules in strong laser fields. In the calcula- 

Sec. HI. A summary and conclusions are given in Sec. IV. 

tions presented here nuclei are treated as classical point 
particles. The Newton equations for the nuclei are solved 

2 Computational method 

for the force which is calculated using the Ehrenfest's 

theorem |23l24j . Electrons are described by the quantum The Coulomb explosion of the carbonyl sulphide is stud- 
TDDFT. They are exposed to the Coulomb field of clas- icd by TDDFT method. Calculations are performed using 
sical nuclei and the laser field. This approach produces the Octopus code where a real-space, pseudopoten- 
errors when the molecule dissociates under a strong laser tial based treatment of the TDDFT has been developed. 
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Physical quantities are expanded using a mesh in the real to calculate the absorption spectrum as 17(0;) ^ 

space. In addition, the ground state properties of OCS are where d{uj) is the Fourier transform of the induced dipole 

also studied by the pseudopotential and planewave based moment d{t) = j zn{r,t)(Pr, and n{r,t) — |?/;(r,t)|^ is 

Abinit code f28^. The TrouUier-Martins pseudopotentials an electron density. The absorption spectrum of OCS, 

[29] are used in both simulations, and the OCS molecule calculated using the Octopus code, shows that the first 

is positioned in the parallelepiped box with the side of 30 peak is at 7 eV. This is in a good agreement with ex- 

a.u. A Cartesian grid with the constant mesh spacing of periments |33I34| . Simulations of the laser-molecule inter- 

0.2 A is applied in simulations based on the Octopus code, action are done at the nonresonant frequency of 2.5 eV, 

The ground-state of OCS obtained using the density func- i.e., below the first absorption peak. Experimental studies 

tional theory is perturbed by the laser pulse. The time- of the Coulomb explosion in OCS have been performed 

dependent Kohn-Sham equations are integrated using the at the nonresonant frequency of A = 790 nm = 1.57 

time step of At = 0.002 h/eV — 1.32 x 10"^^ s. The time eV |13I14| . The length gauge, pulses with a cosine enve- 

evolution of up to 52.8 fs is followed. This is a long time lope, a total length of 6.6 fs, and maximal intensities of 

for the TDDFT method. The approximated enforced time- h ^ 1.3272 x lO^^ w/cm^ and I2 = 1.3272 x lO^^ W/cm^ 

reversal symmetry algorithm is applied to treat the evolu- are applied in these simulations, 
tion operator in the time-dependent Kohn-Sham problem 

|30| . The adiabatic local density approximation exchange 3 Rcsults and disCUSSion 

functional in the parametrization of Perdew and Wang is 
chosen [31 . In this work the velocity Verlet algorithm is 
used to follow the classical time evolution of atomic cores. 
The OCS molecule is oriented along the z axis and exposed 
to the linearly polarized laser light along the x, y, or z axis. 
Both 0-C-S and S-C-0 orientations towards an incoming 
laser pulse are studied for the light polarized along the 
molecule axis. The absorbing conditions are used at the 
boundaries of the box S2] . 



The ground state density functional theory calculation of 
the OCS molecule produces a good agreement with exper- 
imental results. The molecule is linear, the S-C distance 
is 1.548 A, whereas the 0-C one is 1.148 A. The mea- 
sured values of these distances are 1.561 A and 1.156 A, 
respectively [35]. The results for changes of the geometry 
when the laser field is polarized perpendicularly to the 
initial molecule axis are shown in Table I. The C— S dis- 
tance starts to change at 3.96 fs, the O— C distance at 
Electronic excitation spectra of carbonyl sulphide have 9.24 fs, and the angle Z (O— C— S) at 14.52 fs, for the 
been measured from 5 to 360 eV using dipole spectroscopy laser intensity Ii. The O— C distance first increases up to 
[33134] . The first absorption peak is rather broad and 1.209 A (at 19.14 fs), and then decreases. For example, the 
its maximum is at 7-8 eV. Using TDDFT it is possible O— C distance is 1.164 A at 26.40 fs. The C— S distance 
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increases up to 1.689 A at 31.68 fs and then decreases. 



The results for the intensity h = 1.3272 x 10^^ W/cm^ 
also show that the distances do not change monotonicaUy 



1 

1.6 



1.2 - 



in time. The distances between atoms when the laser field u , . 

S 

is polarized along the initial molecule axis are presented 
in Table 11. There the molecule is linear. In this Table 
the results for two orientations of the molecule towards a 
laser pulse are shown. It is calculated that data for these 
two orientations differ only quantitatively. The results for 
the S-C-0 orientation are presented in examples for an 
interaction with the z-polarized laser light in Figs. 2-5. 
Figure 1 shows an example of the change of the O— C and 




20 30 
Time (fs) 

Fig. 1. The time evolution of atomic distances after an in- 
teraction with the laser light of h = 1.3272 x 10^^ W/cm^ x 
polarization. 



C— S atomic distances during the 52.8 fs of the time evo- , , , . . .,, . . ,, 

elongated m comparison with equilibrium ones m the pro- 

lution. Similar nonmonotonic changes of the internuclear „ , , i ; i ^ i i 

cesses ol an enhanced multiple ionization and the Coulomb 

explosion [57] . 



distance have been recently measured and calculated in 
the molecule study by ultrashort laser fields on a time 
scale of 3 ps [55] . 



The OCS molecule ejects electrons after an interaction 
with a strong enough laser pulse. Figure 2 shows how the 

In general, the angle between bonds changes when the electronic charge of this 16- valence electron molecule (as 

laser polarization is along directions which are perpendic- represented by the pseudopotential model) changes when 

ular to the molecular axis. OCS preserves its linear geom- OCS is irradiated by a laser pulse. As shown in Figs. 2(b), 

etry when the laser light is polarized along the molecule charge states of up to 4-3 are observed under the field of 

axis. However, bending of OCS induced by pulses polar- ^ 10^^ W/cm^. It is also found that for all laser light 

ized along the internuclear axis may occur on a longer time polarizations the electronic charge in the box does not 

scale, i.e., after 52.8 fs. The observed maximum of the an- change (within the simulation time of 52.8 fs) under the 

gle between bonds is 23.6° and it is calculated for the laser field of Ii = 1.3272 x lO^'^ W/cm^. The results show 

x-polarized pulse of intensity I2 = 1.3272 x 10^^ W/cm^. that the z-polarized laser light is more efficient in removal 

As in experiments [13114] . it is found that the O— C bond of electrons from the OCS molecule. In this simulation 

is stronger than the C— S bond, and that both bonds often only the overall electronic charge of the molecule is cal- 

strongly increase after an interaction with a laser pulse, culated by integrating the electron density in the simu- 

It is already known that atomic bonds arc substantially lation box. However, in the experimental work and cor- 
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Table 1. The time evolution of the structure for laser intensities h = 1.3272 x 10^^ W/cm^ and h = 1.3272 x 10^^ W/cm^, 
using the polarization direction perpendicular to the initial axis of the molecule (x polarization). The distances d (here, in Table 
II and Fig. 1) are geometrical separations between atoms (including the cases where they are not chemically bonded). 





h 






h 






time [fs] 


d(O-C) 


d(C-S) 


/ (O-C-S) 


d(O-C) 


d(C-S) 


/ (O-C-S) 


2.64 


1.148 


1.548 


180° 


1.148 


1.548 


180° 


7.92 


1.148 


1.560 


180° 


1.152 


1.572 


180° 


13.20 


1.172 


1.576 


180° 


1.184 


1.616 


179.2° 


18.48 


1.204 


1.604 


178.7° 


1.207 


1.688 


170.8° 


26.40 


1.164 


1.668 


166.9° 


1.088 


1.824 


171.7° 


31.68 


1.049 


1.689 


171.7° 


1.012 


1.865 


176.9° 



Table 2. The time evolution of the structure for the laser intensity I2 — 1.3272 x 10^^ W/crn^ and a laser pulse polarized 
along the internuclear axis (z polarization). The data for two orientations (S-C-O and O-C-S) of the carbonyl sulphide molecule 
in the laser field are presented. 
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O-C-S 






time [fs] 


d(O-C) 


d(C-S) 


Z (O-C-S) 


d(O-C) 


d(C-S) 


Z (O-C-S) 


2.64 


1.148 


1.548 


180° 


1.148 


1.552 


180° 


7.92 


1.184 


1.548 


180° 


1.192 


1.568 


180° 


13.20 


1.260 


1.572 


180° 


1.272 


1.600 


180° 


18.48 


1.336 


1.648 


180° 


1.288 


1.660 


180° 


26.40 


1.288 


1.844 


180° 


1.260 


1.784 


180° 


31.68 


1.176 


2.008 


180° 


1.252 


1.880 


180° 



responding Monte Carlo simulations several {Qi,Q2,Q3) 
fragmentation channels have been proposed, where the la- 
bels are as in (0^1 + , C"3=+, 5*^^'+) [T3fT4] . For example, it 
has been found that (3, 3, 5) is a channel with the highest 
ionization under the field of 2± 10^^ W/cm^. 

Several typical examples of the kinetic and potential 
energy are presented in Figs. 3 and 4. It is well known that 
positively charged ions form in the process of the laser- 



molecule interaction. The Coulomb repulsion of the ions 
and an initial binding energy of the molecule are stored 
as the potential energy. The repulsive Coulomb potential 
energy is strong and it is rapidly released as a kinetic 
energy of the molecule fragments. The results show that 
energies change after a time delay of at least 2 fs. Regu- 
larly positioned peaks are clearly visible in Fig. 3(b). The 
peaks are less pronounced under conditions of Fig. 3(a,c). 
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(a) 



10 20 30 40 
Time (fs) 



starts to vibrate and the result is an oscillatory exchange 
of kinetic and potential energy. The potential energy (Fig. 
4(c)) is much larger in an absolute value, and therefore the 
peaks are more visible in the kinetic energy (Fig. 3(b)). It 
has been measured for OCS perturbed by the laser light 
that the kinetic energy release varies with the ionization 
channels [14 . This study has been performed using laser 
pulses focused to ~ 10^^ W/cm^ and the kinetic energies 
between 4.1 eV and 105 eV have been measured. In the 
ground state of the molecule its potential energy is neg- 
ative. When the molecule is perturbed by the laser light, 
the potential energy increases. Figure 4 shows that the 
potential energy is always negative. Therefore, the system 
preserves at least the part of its initial binding energy. The 
time evolution of the electric dipole moment dz{t) of the 
molecule is calculated and presented in Fig. 5. When the 
light is polarized along the molecule axis, the change of an 
electric dipole distribution in the first several fs is espe- 
cially intensive and approximately follows the oscillatory 
Fig. 3. The time evolution of the kinetic energy of OCS after features of the laser pulse. 




(b) 



10 20 30 40 
Time (fs) 



Fig. 2. The time evolution of the electronic charge after an 
interaction with the laser light: (a) h = 1.3272 x 10^'' W/cm^, 
X polarization, (b) I2 = 1.3272 x 10^^ W/cm^, z polarization. 



an interaction with a laser pulse: (a) /i — 1.3272x10^^ W/cm^, 
z polarization, (b) I2 = 1.3272 x 10^^ W/cm^, x polarization, 
(c) h = 1.3272 x 10^^ W/cm^, z polarization. 



The electron localization function (ELF) has become 
an important tool for a description of chemical bonds in 
molecules and solids |38I39| . The results of many studies 
Therefore, the peak shapes and positions are intensity and have been shown that, in comparison with electron densi- 
polarization dependent. These peaks in the time evolution ties, the ELFs better describe the properties of localized 
of the kinetic energy and similar ones for the potential en- and delocalized electrons in atoms and their aggregates, 
ergy (Fig. 4(a,b)) show that an intensive redistribution of Recently the time-dependent electron localization func- 
the energy occurs in the OCS molecule after an interac- tion suitable for an analysis of the creation and breaking 
tion with a laser pulse. Under conditions of Fig. 3(b) the of bonds has been introduced [H]. The three-dimensional 
laser field caused the molecule to bend. Then, the molecule representations of the time-dependent ELF isosurfaccs for 
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Fig. 6. The time-dependent electron localization function ear morphology after 52.8 fs of the time evolution when 
(ELF). An initial ground-state structure is shown in (a). The it is oriented along the laser polarization direction. An 
same molecule after an interaction with the x-polarized laser intensive redistribution of charges and energies, and the 
pulse of the intensity h = 1.3272 x 10^^ W/cm^; (b) 0.66 fs, stretching of bonds are observed. These results are in a 
(c) 2.64 fs. good agreement with experimental studies of the carbonyl 

sulphide molecule in a strong laser field [13|14| . However, 



OCS are presented in Fig. 6. The comparison of ELFs in 
the ground-state (Fig. 6(a)) and after 0.66 fs (Fig. 6(b)) 
shows that the instability of the OCS molecule first oc- 
curs at the sulfur atom (the upper left part of the ELF 
shown in Fig. 6(b)). Figure 6(b) represents the point of 
the system time evolution when the distances are d(0- 
C)= 1.148 A, d(C-S)= 1.548 A. Further analysis of the 



experiments have been carried out on a much longer time 
scale. TDDFT is an ab initio computational method which 
is able to simulate processes in materials at nanoscales in- 
duced by short and strong laser pulses. 



This work has been carried out under the HR-MZOS project 
"Electronic Properties of Surfaces and Nanostructures" . I am 
time-dependent ELF confirms that the electron ejection ^^^^^^^^ ^j^^ University Computing Center SRCE for their 



begins around the S atom. Figure 6(c) shows several blobs j^gjp g^^^j computer time 
of a negative charge. The corresponding geometry of the 
molecule shown in Fig. 6(c) is d(0-C)= 1.148 A, d(C-S)= 
1.560 A. After 2.64 fs electrons are already ejected from 
the molecule. 
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